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Abstract. In this work we introduce a deterministic scheme of synchronization of linear and nonlinear 
cellular automata (CA) with complex behavior, connected through a master-slave coupling. By using a 
definition of Boolean derivative, we use the linear approximation of the automata to determine a function 
of coupling that promotes synchronization without perturbing all the sites of the slave system. 

1 Introduction tems has been addressed. 

There is a remarkably small number are cases where the 
The spontaneous and synchronous evolution of extended model is represented by a totally discrete system such as 
systems or parts of them is an increasingly frecuent phe- CAs |16ll7llg|19l21l22ll9t^ . The before mentioned refer- 
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, nomenon mdisciplmesrangmg from biology to meteorology [T^aililiPg] far as we know, the only cases. Particularly, 
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In several of these cases, the systems are made up of fi- in [16] the authors study the synchronization between two 

nite, discrete and interrelated parts. Nevertheless, the ma- stochastically coupled CA with complex behavior, where 

jor models used to study the synchronization of extended the coupling mechanism is characterized by a probability 

systems involve continuous systems represented by Partial p_ in [17 , the authors study the synchronization between 

Differential Equations [l0|ll|12|13|il] or partially discrete randomly coupled totalistic (hnear) CA; as in the previ- 

models given by Coupled Maps Networks, see for instance qus reference the coupling is characterized by a probability 

[H] and references therein. In both cases the problem of p and the synchronization threshold Pc is related to the 
synchronization of identical and non-identical chaotic sys- 
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maximum Lyapunov exponent of the automata, defined 2 Linear approximation of nonlinear cellular 

in analogy witti continuous dynamical systems. Along the 
same line is the work of Sanchez et. al.[TH]- In this work the 
same stochastic coupling scheme is used to search for sta- 
ble patterns in the evolution of the automata. Rouquier in 
|19| carries out an exhaustive experimental study of syn- 
chronization in CAs using the coupling function proposed 
by [12]. Finally, in a very recent publication[5D], F. Bag- 
noli et. al. use the Jacobian matrix of the CA to estimate a 
probability that later they use to design a coupling scheme 
to reach synchronization or achieve control. 
On the other hand, in references |21I22) the authors de- 
rive conditions for synchronization using a deterministic 
coupling, but only in the case of totalistic automata. 



automata and synchronization 

Cellular automatas are dynamical systems where space, 
time and state variables belong to discrete sets. More for- 
mally, a CA is a 4-tuple (£,5, A/", /), where £ is a d- 
dimensional lattice of cells, S denotes the set of the k 
accessible estates of an individual cell. A/" is a mapping 
defining the neighborhood of each cell and / is a local 
transition rule, that gives the state of the «-th cell in the 
time n + 1 given state of it neighborhood in the time n, it 
is 

4+1 = /(A(^)) (1) 



The aim of this work is twofold. Firstly we introduce a de- 
terministic scheme for coupling Class III, linear and non- 
linear CAs, coupled by a master-slave scheme. The strat- 
egy is based on a sort of linear approximation of the au- 
tomata by means of a definition of Boolean derivative _23._. 
And secondly, we propose criteria to minimize the inten- 
sity of the coupling, in the sense of the number of coupling 
sites. 

The article is organized as follows: in section 2 we intro- 
duce briefly the CA, the notion of the Boolean derivative, 
set the problem of synchronization between CAs and es- 
tablish conditions for synchronization in unidirectionally 
coupled CAs. In section 3, we show the performance of the 
strategy using examples of Class III nonlinear CAs, and 
finally in section 4 we make some concluding remarks. 



Despite it simplicity, CAs exhibit an unsuspected abil- 
ity to emulate, from qualitative 24 and quantitative [25] 
points of view, complex extended systems and offer an 
alternative to model a great variety of systems, such as: 
fluid dynamics, forest fires, ecosystems or epidemics. 

Here, we shall restrict ourselves to the case of one- 
dimensional binary CAs of N cells with neighborhood of 
size (r) equal to 3, defined as elementary by Wolfram in 
[26| . In these models, £ is a linear array where each cell 
adopts one of two possible states in S = {0, 1} and the 
dynamical rule is the form s^^^ — /(■^n"^; •Sn: ■s^+^), it 
means that, A/'(i) {sjj"^, sjj, s^+^}. In this case, there 
are k'' = 256 possible rules /. If {bi} is a binary string, 
corresponding to the result of applying / to each different 
set of 3-tuples (s'~^, s', s'+^) sorted in descending order, 
it is possible to use the Wolfram's naming scheme [2S] to 
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assign, to each rule, a label given by the decimal represen- where © and denotes the binary operations OR exclu- 
tation of the binary sequence {bi}. sive and AND respectively and K is a coupling matrix 

Rules of elementary automatas were classified qualita- that define the coupling, this is, the number and location 
tively by Wolfram[26] according to its asymptotic behavior where the master system perturb to the slave system, 
as Class I to IV, if the automata evolution asymptotically The evolution of the difference automata defined as 
goes to: a constant state, an isolated periodic segments, An+i = Sn+i (BSn+i, is given by: 
chaotic regime or isolated chaotic segments, respectively. 

Zi„+i = F(S„) ® F(S„) © K (S„ © S„) 

In contrast with the continuous case, the elements of 

« F{Sn) © F(S„) © J(S„) /!„ © K Z\„ 

the Class III can be linear or nonlinear CA. The linear are 

the rules that can be written as ^ (^(S") ® ^) ® ^" 

„, ,_i , ,+is \ - „ , where i^(S„) is linearly approximated as in Frf and the 

/« - mo(i(2) aj e{0,l}. (2) ^"^ ^ - 

Jacobian matrix of F, J(S„) was defined in [53], as 

In the case of elementary CA this class contain 16 elements 

, |^ = F(s^...,s^...)ffiF(s^...,s^■©l,...) (6) 

and the rest are nonlinear. as^ 

For an array of size iV it is possible to define the au- matrix elements- 

tomata state as S„ = {s^j, . . . , s^} whose evolution is 



governed by ds^ 



J.^ - id (7) 



S„+i = i^(S„) (3) 

It is worth noting that this Jacobian is a tri-diagonal 

, „ . 1 , , , • 1 1 , , , , 1 ,. matrix except for two off-diagonal elements in the corners 

where is a global rule induced by the local rule /. 

for periodic boundary condition. This structure could be 

In our problem, we deal with two identical Class III au- 

~ colloquially explained as follows: each element in the prin- 

tomatas S„+i = _F(S„) (the driver system) and S„+i = 

~ . . . . . cipal diagonal represent a cell, an 1 in this diagonal indi- 

F{Sn) (the driven system) with periodic boundary condi- 

cate that this site is affected in the future by a change of 

tions and different initial conditions So and So- 

it in the present; ones in the column indicate that this site 

In analogy with an usual strategy in the continuous 

is affected in the future by actual changes in its neighbors 

case, we connect the systems by mean of an unidirectional 

and ones in the file indicate that changes in this site affect, 

linear coupling that penalize the separation between the 

in the future, its neighbors. 

estates of the driver and driven systems in the form 

From ([5]), it is clear that if K = J(S„) the difference 
^"+1 ~ -^(^n) automata converge to zero and the master-slave system 

S„+i = F(Sn) © K (S„ © S„) (4) synchronizes. Nevertheless, this coupling represents a con- 
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nection between all elements of the master and slave sys- tion of the perturbation over the slave system can drive 

tems by mean of a coupling function that varies in time. its state close to the state of the master system using a 

In order to design a more efficient coupling function (in suitable value of q. 
the sense of the complexity and intensity of this function) 

we suppress the time variation approximating J(Sn) by ^ Numerical Rssults 
J(Sn,q), as: 



In order to show the performance of the strategy we carry 



3{Sn,q)-H{q ^2^J(S„)) (8) out experiments with all Wolfram's Class III CAs with 

where H is the Heaviside function and < g < 1 is a non-constant Jacobian matrix and present as examples 

threshold parameter that represents a kind of tolerance ^^ed in [T^ i.e. automatas: 18, 22, 30. 41, 45, 54, 

in the approximation of the local Jacobian by the aver- ^O^, 110, 122, 126 and 146, with N = 100. We calculate 

age Jacobian. In the particular case of linear automata, normalized Hamming distance H{q) and the fraction 

J coincides with the Exact Jacobian independently of the couphng sites, 

value of q. , w N number of coupling sites 

^ ^c{q) = , (11) 

In this way, the resultant Jacobian could yet have a 

as a function of the parameter q, averaged for 500 ran- 

structure so that its application in ([S]) results in the per- 

domly chosen initial conditions. Here we can easily iden- 

turbation of all sites in the slave system. In order to mini- 

tify three types of behavior, that are shown in the Figures 

mize the number of coupling sites, we will chose K, select- 

(HJ , © and Q : the synchronization occur in a interval of 

ing the coupling sites as the sites where the perturbations 

values of q approximately between 0.5 and 0.7, the syn- 

are more easily spread, and at the same time, are sites 

chronization occurs in an interval of values of q approxi- 

more sensitive to be perturbed by changes in its neigh- 
mate between y 0.0 and 0.75, the synchronization occurs 

borhood. These are sites (c) so that the Jacobian matrix 

in an interval of values of q approximately between 0.0 and 

has 2 or 3 ones in the corresponding row (Jc,i) or column 

0.5. In all cases there is a transition from a high number of 

iJ^,c)■ 

{J ( \ ' h --^ A coupling sites to uncoupling state around the upper limit 

(9) . . . . ^ 

si < 4 '^^ previously mentioned intervals. Figure (14]) shows 

where h is minimum number of coupling sites required to achieve 

Ci \ / 1 \ synchronization averaged over the 500 randomly chosen 

J2 J^+nAq)] E J.a+niq)) (10) 
J / initial conditions. As you can see in this figure, there are 

If F is nonlinear S„+i 7^ J(g)S„ and of course Sn+i 7^ CAs that should be perturbed in almost all sites in order 

3{q)Sm nevertheless as we will see, a reiterative applica- to reach a synchronized state. This is because the matrix 
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K, which determines the number of coupling sites, is the ^ilX^c(g) 
result of averaging 500 realizations started from randomly 
chosen states and it is quite likely to depend on the dif- 
ference between the initial conditions, as suggested in[14j. 




H(q). N,(q) 




Fig. 3. Normalized Hamming distance H{q) (continuous 
line) between asymptotic states of the master and the slave 
systems and Normalized number of coupling sites Nc{q) 
(dashed line) as a function of q, averaged for the case of 
the automatas 45, 54, 106. 



Fig. 1. Normalized Hamming distance H{q) (continuous 
line) between asymptotic states of the master and the slave 
systems and Normalized number of coupling sites Nc{q) 
(dashed line) as a function of q, averaged for the case of 
the automatas 18, 110,126, 146. 
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18 110 126 146 22 30 41 122 45 



Fig. 4. Minimum number of sites necessary to reach syn- 
chronization for all CAs, averaged over 500 realizations. 

As an specific example, we use the automata rule 22: 



n ' n ' n 



^©<0<+^ (12) 



Fig. 2. Normalized Hamming distance H{q) (continuous 
line) between asymptotic states of the master and the slave ^^^^ pjg^^^ gj, ^ evolution of the master-slave 



system when the parameter q — 0.802. In this particular 



systems and Normalized number of coupling sites Nc{q) 
(dashed line) as a function of q, averaged for the case of ^^^^^ numhev of couphng sites was N, = 0.66 
the automatas 22, 30, 41, 122. 



Finally, it is important to note that during the numeri- 
cal experiments we observed that when used as coupling a 
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Fig. 5. Evolution of master, slave and difference automata 
for the rule 22 and q = 0.802. The horizontal axis represent 
the space and the time goes from top to bottom 

constant matrix randomly generated, the synchronization 
is achieved in a small number of attemps. 



4 Final Remarks. 

We present a deterministic master-slave coupling scheme 
to achieve synchronization in linear and nonlinear CAs. 
Here the master, in general, need not act on the the full 
extent of the slave system, in contrast to probabilistic 
schemes where it is necessary to monitor the entire slave 
system all the time, as well as the posibility of disturbing 
all the sites. 

The scheme is valid in linear and non-linear cases. Par- 
ticularly in the case of linear CA, the Jacobian matrix 
coincides with the global rule F ([3]); hence the synchro- 
nization is reached in one iteration using K as defined in 



The strategy presented seems to be useful in the design 
of strategies to control complex extended systems, that 
can be modeled by cellular automatas, such as epidemics 
or forest fires, by a efficient pinning. 
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